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Properties of a scalar a-meson are investigated in the two-flavor Nambu-Jona- 
Lasinio model with the Polyakov loop. Model analysis of the phase diagram of strong 
interacting matter is performed. The temperature dependence of the a — ► tttt decay 
width is studied at the zero chemical potential and near the critical end-point. The 
calculated strong coupling constant g a -mr and the decay width are compared with avail- 
able experimental data and other model results. Nonthermal enhancement of the total 
decay width is noted for the a meson near the critical end-point when the condition 
m a > 2m T is broken. 
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Introduction 

The models of Nambu-Jona-Lasinio type [TJ [2j [3l S] have a long history and 
are used to describe the dynamics and thermodynamics of light mesons. This 
type of models gives a simple and practical example of the basic mechanism of 
spontaneous breaking of chiral symmetry and key features of QCD at finite tem- 
perature and chemical potential [SJ |H1 [5J IH] • The behavior of a QCD system 
is governed by the symmetry properties of the Lagrangian, namely, the global 
SUL(Nf) x SUn(Nf) symmetry which is spontaneously broken to SUv(Nf) and 
the exact SU C (N C ) local color symmetry. On the other hand, in a non-Abelian 
pure gauge theory, the Polyakov loop serves as an order parameter of a tran- 
sition from the low temperature confined phase (Zn c symmetric) to the high 
temperature deconfined phase characterized by the spontaneously breaken Zjq c 
symmetry (PNJL model). In the PNJL model, quarks are coupled simultane- 
ously to the chiral condensate and to the Polyakov loop, and the model includes 
the features of both the chiral and Z/v c symmetry breaking. The model repro- 
duces rather successfully lattice data on QCD thermodynamics. The use of the 
PNJL model is therefore reasonable for investigating the in-medium properties 
of mesons and their decays [TQl [IT] . 

The aim of this work is the investigation of the meson properties and a 
decay near the critical end-point (CEP). In this letter, we discuss the decay 
process a — > tttt at the finite temperature T and chemical potential fi in the 
framework the Nambu-Jona-Lasinio model with the Polykov-loop (PNJL) that 
is believed to describe well the chiral properties and simulates a deconfinement 
transition. Our motivation here is to elaborate these features in a large region 
of the temperature T and quark chemical potential /i, where a non-thermal 
enhancement of pions due to the a — > tttt decay may take place. 
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1 The model and the phase diagram 



We use the two-flavor PNJL model with the following Lagrangian [12j [T3l [T^ 



£p N JL = q^D» -m )q + G[(qqf + (qi l5 rq) 2 \ - It (<f>[A], $[A}; T) , (1) 

where the covariant gauge derivative = — iA^ with A^ — 5^A a , A — 
—iA± (the Polyakov calibration). The strong coupling constant is absorbed in 
the definition of A^. At the zero temperature the Polyakov loop field $ and 
the quark field are decoupled. Here, the quark field q — (u,d), current masses 
77i = diag(m„, mj), Pauli matrices r — a/2 act in the two color flavor space and 
G is the coupling constant. 

The gauge sector of the Lagrangian density ([T]) is described by an effective 
potential U ($[A],$[A];T) = U ($, T) 



W($,$;T) 

rp4 



b 2 (T) 



(2) 



where 



b 2 (T) = a Q + ai 



2 



0.3 



(3) 



The parameter set is obtained by fitting the lattice results in the pure SU(3) 
gauge theory at T = 0.27 GeV OH E] and is given in Table [U 



a 


ai 


a2 


a3 


b 3 


fo 4 


6.75 


-1.95 


2.625 


-7.44 


0.75 


7.5 



Table 1: The parameter set of the effective potential $; T). 



Before discussing the meson properties, one should indroduce the gap equa- 
tion for constituent quark mass should be introduced. For describing the system 
properties at the finite temperature and density the grand canonical potential 
in the Hartree approximation is considered []2J [H] 



n(M,m,T,p) = U ($, #; T) + N f (m ^ - 2N C N S 



d 3 p 



(27T) 



d 3 p 



[in N+(E P )+ In N^(E P )] , 



(4) 



where E p is the quark energy, E p = \/p 2 + m 2 , E p = E p =p /z, and 



1 + 3 ($ + $e _/3i ^ 
1 + 3 ( $ + <£>e~ pE p ) < 



-3PE- 



(5) 
(6) 
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Integrals in Eq. Q contain the three-momentum cutoff A. 

From the grand canonical potential f2 the equations of motion can be ob- 
tained 

and the gap equation for the constituent quark mass can be written as follows: 

d 3 p m 



rn = rn - N f G <qq>= m + 8GN c N f . , 
where f£, are the modified Fermi functions 



(8) 



/+ = (($ + 2$e-^ + )e-^ + + e~^ E+ )N+, 



/»=((* + 2*e-^ )e 



-/3_E _ \„-/3B" 



<I> • 



(9) 



with i?^ = _E =p /i. The regularization parameter A, the quark current mass 
mo, the coupling strength G and physics quantities to fix these parameters are 
presented in Table [2] 



m [MeV] 


A [GeV] 


G [GeV]- 2 


Fa [GeV] 


mir [GeV] 


5.5 


0.639 


5.227 


0.092 


0.139 



Table 2: The model parameters and quantities used for their tuning 



The a and n meson masses are the solutions of the equation 

l-2GII Ws (fc 2 )=0, (10) 

where k 2 = m 2 , and k 2 = m 2 in pseudo scalar and scalar sectors, and H ps / S are 
standard mesonic correlation functions [7] 

— 4jj Tr [i l5 r a S(p + k)i l5 r b S(p)] , (11) 

iUAk ' 2) = Jw? Tr [lS[p + k)tS(p)] • (12) 

Both pion-quark gcr W7r (T, /i) and sigma-quark go-^ coupling strengths can be 
obtained from n ps / s : 

-a \ dll^fc 2 ) 1 /,ox 

As is seen from the gap equation ([5]), the quark condensate < qq > de- 
fines completely the quark mass in a hot and dense matter. This correlation is 
clearly demonstrated in Fig. [TJ where the temperature dependence of the order 
parameters of the chiral condensate and the Polyakov loop, as well as the sigma 
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Figure 1: The temperature dependence of the chiral condensate and Polyakov 
loop (left panel) and particle masses (right panel) at /i = within the PNJL 
model. 




Figure 2: The phase diagram within the PNJL model. The solid lines denote 
the first-order transition boundary, the dotted line is the second-order transition 
and the dashed line is a crossover. The lines with the CEP and TCP points are 
calculated for a finite mass and in the chiral limit m q =0, respectively. 



and pion masses are shown for /i = 0. The pion hardly changes its mass and 
starts to become heavier only near T c ~ 200 MeV, while the sigma mass m a {T) 
decreases as the chiral symmetry gets restored, and eventually m ff (T) becomes 
larger than double quark mass m q at the temperature T\iott ~ 190 MeV. 

Since it is still difficult to extract certain information from the lattice simu- 
lations with the nonzero baryon density, we need QCD models for investigating 
the phase transitions at the finite baryon density. The calculated phase dia- 
gram of the physical states of matter within the PNJL model is given in Fig. [5J 
As is seen, in the real world with nonzero pion mass we have the first-order 
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phase transition at a moderate temperature and a large baryon chemical po- 
tential fiB = 3/i that, with increasing T, terminates at the critical end-point 
{Tcep, Hcep) where the second-order phase transition occurs. At a higher 
temperature T > Tctp we have a smooth crossover. In the chiral limit with 
massless pions there is a tricritical point that separates the second-order phase 
transition at high temperature T and the first-order transition at lower T and 
high /i. For the model parameters chosen (see Table[2]) we obtain Tcep —0.095 
GeV and Hcep —0.32 GeV (cf. with [T5], where the critical temperature and 
chemical potential are calculated with the same parameters). 



2 Decay a — > irir 

It is in order here to mention the significance of the scalar meson a (a chiral 
partner of the pion) in QCD. A model-independent consequence of dynamic 
breaking of chiral symmetry is the existence of the pion and its chiral partner 
cr-meson: The former is the phase fluctuation of the order parameter qq while 
the latter is the amplitude fluctuation of qq. During the expansion of the system, 
the in-medium a mass increases toward its vacuum value and eventually exceeds 
the 2m n threshold. As the a — » tttt coupling is large, the decay proceeds rapidly. 
Since this process occurs after freeze-out, the pions generated by it do not have 
a chance to thermalize. Thus, one may expect that the resulting pion spectrum 
should have a nonthermal enhancement at low transverse momentum. 

To the lowest order in a \/N c expansion, the diagram for the process a — > tttt 
is shown in Fig. [3] The amplitude of the triangle vertex a — > tttt can be obtained 
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Figure 3: The Feynman diagram of the a — > itit decay 



analytically as 

j^L Tr{S(q) r OT S(q + P) S{q)}, (14) 

where IV = 172^" is the pion vertex function and S(q) = q + ra/q 2 — m 2 is the 
quark propagator, a trace is beeng taken over color, flavor and spinor indices. 
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Figure 4: The temperature dependence of the total decay width and the 

coupling strength g a7r7r at the zero chemical potential. The solid and dash-dotted 
lines are the results for the decay width (fT6|) with and without the Bose-Einstein 
factor Ftttt, respectively. 
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After tracing and evaluation of the Matsubara sum one obtains [THl 02] 



(2^) 3 2E q 

o) + m,2/2 - 2m 2 /? 2 

(15) 



2mN c Nf 

(q ■ p) 2 - (2m 2 + 4m 2 )(q ■ p) + m|/2 - 2m 2 E g 2 



(m 2 -4£2)(( m 2 -2(q.p)) 2 -m 2 £;2) 

where are the modified Fermi functions ([9]). The coupling strength 

9<nrv(T, /x) = 2g a g 2 A CTW7r (T, /x), where </„ and are coupling constants defined 
from Eq. i|I5]l 

The decay width is defined by the cut of the Feynman diagram in Fig|3] 
treating the sigma meson as a quark- antiquark system 



r _ 3 flg ffT / 4m 2 

r ™~2 le^y 1_ ^r (16) 

The scalar meson a can decay either into two neutral or two charged pions. 
All these channels are taken into account. The factor 3/2 in Eq.(|16p takes into 
account the isospin conservation. 

For the existing decay a — ► -kit the kinematic factor yjl — Am^/m 2 in (|16|) 
leads to the constraint 7710- < 2m,,. which is T- and /^-dependent. One may 
expect that in the temperature region where this condition is broken the values 
(7 CT7r7r and r^j — ^7T7r will drop to zero (Fig. 0|). In the case /x = considered the 
kinematic condition is broken at the a — > irn dissociation temperature T[[ 
210 MeV. 

The coupling constant ga-mv is about 2.1 GeV (note the scaling factor 1/10 
in Fig|4j in vacuum and stays almost constant up to T ^ 0.22 GeV (at /i = 0) 
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Figure 5: Temperature dependence of a (left panel) and 7r (right panel) me- 
son masses near the critical end-point. Solid lines correspond to the chemical 
potential at the critical end-point ficEP- Dot-dashed and dashed lines are the 
PNJL results for ficEP + 10 and p»cep — 10 MeV, respectively. The dotted lines 
correspond to the mixed phase. 




Figure 6: Left panel: The temperature dependence of the coupling strength 
ga-TTTr at three values of the chemical potential /i = fxcEP, hcep-10 MeV and 
Hcep+W MeV. Right panel: The total decay width r <J _ ! . Tr7r at the same values 
of chemical potentials. The total decay width F n7r and the width IV^, 

Eq. (fT6| are plotted by the solid and dashed lines, respectively, for the same 
three values of /x. 



and then it drops to zero at m a — 2m T . The experimental value, extracted from 
the J/V> decays, g aiVK = 2.0l ; g GeV [T5j is in reasonable agreement with our 
result. It is of interest to note that the quark-meson models predict g Gvr: =1.8 
GeV [19] and 1.8^q'| GeV [20], and the linear sigma model gives g a -K-K = 2.54± 
0.01 GeV gT] . 

The total a decay width was measured recently in two experiments [15] [201 
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22J. The Beijing Spectrometer (BES) Collaboration at the Beijing Electron- 
Positron Collider reported evidence of the existence of the a particle in J/ijj 
decays. In the tt + tt~ invariant mass spectrum in the process of the J/ip — ¥ 
aoj — > it + it~oj they found a low mass enhancement, and the detailed analysis 
strongly favors ++ spin parity with a statistical significance for the existence 
of the a particle. The BES measured values of the a mass and total width 
are [HHD] 

m a = 390+^° MeV, = 282+^ MeV. (17) 

The E791 Collaboration at Fermilab reported on evidence of a light and broad 
scalar resonance in the nonleptonic cascade decays of heavy mesons [TH] • It was 
found in the Fermilab experiment that the a meson is rather important in the 
D meson decay D — > Sir generated by the intermediate c-resonance channel 

m a = 478123 ± 17 MeV > T<r-+™ = 324±^ ± 21 MeV. (18) 

These experimental values should be compared with T^^^ —190 MeV at 
T = fi =0 (see Fig[4|. One can additionally take into account the Bose- Einstein 
statistics of final pion states by introducing the factor F n7r = (1 + fsi^)) 2 Q2] 
into the total decay width (|16l) . where the boson distribution function fsix) = 
(e x l T — In contrast to the kinematic factor yl — 4mJ/mJ, this pion 

distribution function tends to increase the width. Near the Mott temperature 
r<T->7r7r ~250 MeV at [i =0. In fact, the numerical calculation shows that 
r (T7r7r (T) decreases as T goes up, and eventually vanishes at a high temperature. 
The measured widths (fl7 jl . (|18p are somewhat higher then that in our model at 
T =0. It is noteworthy that the measured a masses 390±3g [H] and 478±23 [22] 
are noticeably smaller than in our model m a «620 MeV (to be fixed by the 
model parameters) while the total decay width depends strongly on m a , see 
Eq. (|16l) . The quark- meson models give the decay width that does not differ 
essentially from our estimate: To-^^ =173 [19] and 149.9 [23] MeV though the 
used sigma meson mass is close to experimental ones, being m a = 485.5 and 
478 MeV, respectively. 

The decay width was really not studied at the nonzero baryon density, in 
particular, in the region near the critical end-point. For the finite fj, the coupling 
strength and meson masses behave in a nontrivial way. Both a and 7r mesons 
suffer a jump in the region of the first-order phase transition (see curve for /Lt+10 
MeV in Fig[5]) which ends at the critical end-point hcep and then they change 
continuously for < /j-cep, where the crossover-type phase transition occurs. 

According to (fT6| . this mass behavior defines the total decay rate of a a 
meson at // ^ 0, as shown in the right panel of Fig. [6] As can be seen, the 
T-region of the decay enhancement becomes narrower with increasing /i. The 
total widths To--^^ and r ' a ^ wn F wn tend to grow with temperature exhibiting a 
narrow maximum due to particularities of the term in (|16p. The inclusion 
of the F^tt factor enhances this maximum, but this effect becomes smaller when 
one moves by A[i toward larger /i from hcep because the pion mass increases 
above hcep (see Fig[5])- If the chemical potential /i = /j-cep, the temperature 
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at which coupling strength g a ^ or total decay width Y a ^,^ drops down to zero, 
the temperature is just T = Tcep- If M increases (decreases) with respect to 
Vcep by A/i ~10 MeV, the temperature decreases (increases), respectively, by 
about 30 MeV from Tcep- Maximal values of 1^^^ near the critical end-point 
are larger than that for /x =0 by a factor about 3. It means that a a meson lives 
shorter in the dense baryon matter. The shape of g^^iT) is insensitive to /x 
around (J-cep- The full decay width Ya-^^ exhibits some wide maximum near 
the temperature Tcep- This maximum is enhanced due to the factor and 
gets more pronounced for smaller /x. 

Concluding remarks 

The two-flavor PNJL model that reasonably describes quark-meson thermody- 
namics at finite temperature and chemical potential is applied to calculate the 
a — > tttt decay in this medium. The emphasis here is made on the behavior near 
the critical end-point. At fi =0 and T — >0 the coupling constant g a ^ v =2.1 GeV 
and the total a decay width Y a ^ V7r «190 MeV are in reasonable agreement with 
both available experimental data and quark-meson model estimates. At finite 
quark chemical potential near the critical end-point the Y^^^^ width shows 
a sharp maximum coming from a particular behavior of the coupling strength 
g a7r ir- The sigma mesons live here a shorter time than in the baryonless mat- 
ter. The rapid decrease of at a high temperature is due to the phase 
space factor (see Eg . (fl6|) ) . The account for the Bose-Einstein statistics in the 
final state pions (the factor F V7r ) results in the appearance of some nonthermal 
maximum of the decay width near T and /x at which the kinematic condition 
m a < 2m T is broken. This width enhancement is about ~ 20% at iicep an d is 
negligible if one moves to /x — IJ>cep + A/x. 

The presented results are obtained in the first order in the 1 /N c expansion. 
In a more realistic case the a — u> and a — A\ mixing can affect noticeably the 
considered quantities, especially for the /i ^0 case [H]. However, it corresponds 
to higher orders in 1/N C . 

The measurements of nonthermal enhancement of pions might be considered 
as a signature of chiral phase transition. However, it is a difficult experimental 
problem since the a life time is very short, and the pion contribution from the 
resonance decay should be separated carefully. So a more elaborated analysis is 
needed. 
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